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(U ■ Abstract 



Effects of charge-density fluctuations on a phase behavior of the restricted primitive model 
(RPM) are studied within a field-theoretic formalism. We focus on a A-line of continuous transitions 



j^ I between charge-ordered and charge-disordered phases that is observed in several mean-field (MF) 

^ ■ theories, but is absent in simulation results. In our study the RPM is reduced to a (j)^ theory, 

P^ . and a fluctuation contribution to a grand thermodynamic potential is obtained by generalizing the 

o, 

O , Brazovskii approach. We find that in a presence of fluctuations the A-line disappears. Instead, a 

fluctuation-induced first-order transition to a charge-ordered phase appears in the same region of 
i/^ • a phase diagram, where the liquid - ionic-crystal transition is obtained in simulations. Our results 

m : 

en . indicate that the charge-ordered phase should be identified with an ionic crystal. 
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I. INTRODUCTION 



Molten salts, ionic liquids or electrolytes can be described by the restricted primitive 
model (RPMj, where impenetrable hard cores of diameter a carry charges with equal mag- 
nitude e jjj, |2]. In the continuum-space RPM a separation into uniform ion-dilute and 
ion-dense phases with an associated critical point occurs at low densities, a transition to 
an ionic crystal of the CsCl type occurs at intermediate densities, and at high densities the 
fee crystal is stable Q]. The above phase-behavior was confirmed by recent simulations j^. 
n the last decade a very intensive debate was focused on critical properties of the RPM 



less attention 



10 



11 



12 



la 



M 



15|. Phase transitions to crystalline phases drew much 



ITJ I until very recently 



Q 



In addition to the phase transitions found in simulations, a line of continuous phase 



transitions (A-line) was found in theoretical studies |l|, |3|, llSL Il9l. 2G, |2l|, |2^ |23|, except 
from the mean-spherical (MSA) and related approximations yJ,|2J,|25|. Along the A-line a 
decay length of a charge-density correlation function, which exhibits exponentially damped 
oscillations on the length scale ~ a, diverges. In some theories the A-line is_separated from 
a first-order transition by a tricritical point (tcp) 
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y. In Ref. y 



this line was 
just rejected as an unphysical solution. Indeed, a location of the A-line on a phase diagram 
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depends strongly on a regularization of the Coulomb potential inside the hard core 
This fact may indicate that the A-line is an artefact that results from approximations made 
in different theories |lJ, |2i|. On the other hand, in Ref. [iSf it was conjectured that a 
divergent correlation length is a signature of a crystallization. No quantitative arguments 
supporting the above conjecture were given, however. Thus, a role of the A-line in the 



approximate theories 




(all of them of a mean- field (MF) type). 



and its existence in the RPM, remained unclear 

Renewed interest in the whole phase diagram of the RPM, especially in the A-line and the 
tcp, is motivated by recent results obtained for the lattice RPM (LRPM), where positions of 
ions are restricted to sites of different lattices. On lattices with different symmetries, and/or 
with a lattice constant a corresponding to different values of a /a > 1, the ions form different 
periodic patterns at low temperatures T and/or at high densities p. Different patterns 
correspond to different charge-ordered phases. Transitions between the high-temperature, 
charge-disordered phase and the charge-ordered phases are either continuous or first-order. 



depending on details of a lattice structure 
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34 



35, 
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37|, 



38|. In particular, 



on a simple cubic lattice (sc) with a /a = 1, only an order- disorder transition to a phase 
with two oppositely charged sublattices occurs; this transition is continuous for p > ptc, 
where ptc denotes density at the tc£. The phase separation into dilute and dense, uniform 



phases is only metastable [13^ 
by the MSA for the LRPM 




2Ct l27l |28[ . Note that no continuous transition is predicted 
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J39|, in an obvious disagreament with simulations 
37l | and exact theoretical predictions [39|. The two types of the charge-ordered - charge- 
disordered transition are shown in Fig.l. According to recent simulations jj], the transition 
lines between the liquid and the CsCl crystal are very similar to the thick lines shown in 
Fig.l. Note that in contrast to close-packed crystals, the transition density shows significant 
dependence on temperature. The above observations raise a question on a relation between 
the A-line in the continuum space and the charge-ordered - charge-disordered transitions on 
the lattice. 

In this work we study effects of fluctuations on the A-line within the field-theoretic descrip- 
tion developed in Ref . |l^ . On a MP level of this theory the phase diagram for each version 
of the RPM is the same as on the sc lattice with a/a = 1 (thin lines in PigJJ. Namely, only 
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When 



the order- disorder transition that is continuous for p > pt^ is present 

fluctuations are included within the field-theoretic approch initiated by Brazovskii 40|, in 



some lattice systems the order-disorder transition becomes fluctuation-induced first order 
Um (tHic. Un. .„Fi^,. T,e o.de. of the Ua^Mo,. a,.ee. with ...ation .e.u,t. 

, I22I, |2S|- In Ref. [2a, I2JI arguments were given that in the 



for all considered cases [Sll, l2d, Oil, 1^ • In Ref. | 

continuum-space RPM the order- disorder transition becomes fluctuation-induced first order 

as well. 

Except from the order of the considered transition, its location on the phase diagram is 
of major importace for an identification of the charge-ordered phase. On the MP level of 
our theory the order-disorder transition occurs at low densities and high temperatures, and 
the separation into uniform ion-dilute and ion-dense phases is suppresed. Beyond MP, and 
under the assumption that the order-disorder transition is moved away by fluctuations, the 
considered field theory jl3| predicts that for low densities the phase separation into ion-dilute 
and ion-dense phases occurs. The associated critical point belongs to the Ising universality- 



class 
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15|, |26|, i n ag reement with the earlier theoretical arguments by Stell 



with recent theory [14!], experiments 
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and simulations 
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FIG. 1: Schematic representation of the order-disorder transition in the LRPM. Thin dashed- and 
soHd lines represent continuous and first-order transitions respectively that were found on the sc 
lattice with a ja 



1 



m 



m 



31 



m 



37t |. The dashed line is a lattice-analo g of the A- 
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35|. 



Thick solid lines represent the first-order transition that occurs for a ja = \/2, 2 
The transition shown by the thin lines can be continuously transformed to the transition shown 
by the thick lines when additional nearest neighbor repulsion J is present. For small values of J 
the diagram is shown by the thin lines. When J exceeds a certain value, Jq, the dashed line splits 
into two lines that move away when J increases 



M 



, and the first-order transition that occurs 



for large values of J is represented by the thick solid lines. According to simulation results 
the shape of the liquid-CsCl crystal two-phase region in continuum-space RPM is similar to that 
shown by the solid lines. Other transitions that occur in some versions of the LRPM and in the 
continuum space at low- and at high densities are not shown. T and p are in arbitrary units. 

is necessary to verify if the fluctuations may lead to a shift of the phase boundaries of 
the charge-ordered phase from the phase-space region where the gas-hquid separation takes 
place, to the phase-space region where the CsCI crystal is stable, to make the field-theoretic 
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15l l2q complete, and to identify the 



arguments in favor of the Ising universality class 

charge-ordered phase with the CsCl crystal. This is a purpose of our work. 

Our work is based on the Brazovskii theory, which turned out to be sucessful in a de- 
scription of phase transitions and structure of soft-matter systems J47|, U^ Ii9| • Analogous 
theory for hard crystals has not been developed yet. The important common feature of the 
soft- and ionic crystals is that the periodic ordering is not a result of close packing, but fol- 
lows directly from interaction potentials, or effective, state-dependent potentials that favour 
periodic structures for any density. Since the leading physical-mechanism that induces the 



periodic ordering of soft- and ionic crystals is similar, we expect that the Brazovskii approach 
is an appropriate description of ionic crystallization. 

In sec. 2 the field-theoretic description of the RPM is described, our approximations are 
discussed, and notation is fixed. In sec. 3 we derive approximate expressions for the grand 
potential with the fluctuation-contribution included. The following section is devoted to the 
results obtained for the order-disorder transition. The last section contains a short summary 
and a discussion. 



II. FIELD-THEORETIC DESCRIPTION OF THE RPM 



Field theory for the RPM that is considered in this work was derived in Refs. 



u 



26 



35|. In this section we summarize the key steps of the derivation, discuss assumptions 
and approximations, and fix our notation. We consider local deviations from the uniform 
number- and charge densities, ri{-x) = p*(x) — Pq = p+(x) + pl(x) — p^ and 0(x) = p+(x) — 
pl(x) respectively. p^(x) and pl(x) correspond to a local number- density of cations and 
anions respectively, and Pq is the most probable number density of ions. Asterisks indicate 
that all densities are dimensionless, and the unit volume is a^, where a is the core diameter. 
is the charge density in e/a^ units, e is the charge. We focus on systems that are globally 
charge-neutral, 

' '(x) = 0, (1) 



where in this paper we use the notation J = J dx. Deviations from equilibrium, uniform 
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distributions of ionic species are thermally excited with the probability density 

p[<P,p*]=E-'exp{-(3Q''^[<P,p*]), (2) 

where S is a normalization constant, and in our theory f2*^^ is approximated by 
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Q''^[<P, p*] = F,[0, p*] + U[<P] - p /" p(x). (3) 

Jx 

fi is the chemical potential of the ions, F^ = J fh is the hard-core reference-system Helmholtz 
free-energy of the mixture in which the core-diameter a of both components is the same. For 
the continuum RPM we adopt the Carnahan-Starling (CS) form of fh in the local-density 



approximation, 

/?A(P , </>) = ^- log (^- j + ^- log (^^- j - p + p ^^-^, (4) 

where p* = rj + p^, for p* = pg, the f2^^^[0, p*] assumes a minimum, and s = np* /6. Finally, 
the energy in the RPM is given by 

m<p] = f^ / / ^(ix'-x| - i)^w^ = ^ fmnm-^), (5) 

'^ JxJx' |X — X I ^ Jk 

where J^ = J dh/ {27r)^ . Contributions to the electrostatic energy coming from overlaping 
cores are not included in (jSJ. We should note that the regularization of the Coulomb 
potential for r < a is to some extent arbitrary; in particular, in Refs. J2l|, |22, l2i| different 
regularizations were chosen. Here and below x = |x| is measured in a units. (3* = 1/T* = 
j3e'^/{Da) is the inverse temperature in standard reduced units; D is the dielectric constant 
of the solvent. V{k) = Att cos k/k'^ is the Fourier transform of V{x) = 6{x — l)/x, and k is 
in cr~^ units. From the minimum condition for f2*^^[0, p*] we obtain the relation between 
Pq and the intensive parameters, 

si8-9s + 3s^ 
/3/i = logpo + (TZ^3 • (6) 

The fields and rj occur with the probability 0, where the functional f2*^'^[0,pQ + rj] 
consists of a constant term Q^'^^[0,Pq] which is irrelevant, and of the term that depends on 
and ?7, 

An^'^[<t>, r/] = n''^[<f), pl + v]- ^''^[0, Po] = fi2[0, v] + ant[0, v]- (7) 

The boundary of stability of AQ'^'^^[(f),ri] is determined by the Gaussian part, 

1 r dk 



Pn, 



2 J (27r)3 



C" (A;)0(k)0(-k)+7o,2r](k)r/(-k) 



(8) 



where 



and 



^lW=Po +/5*V^W, (9) 



d'^(3fh _ 1 + 4s + 4^2 - 4s3 + s^ 

''dp^\p*=Po - (1 - sYpI 






when the CS reference system is used. The boundary of stabihty, C^^{kb) = 0, occurs along 
the A-hne 

T* = -V{h)pl ^ 1.61p*, (11) 

where fc;, ~ 2.46 corresponds to the minimum of V{k) |23L l26l. l85l| . 

The last term in Eq.((7j) is local, and can be written as 

/3fii„40,r/]= f pujint{(p{x),ri{x.)), (12) 

with 

2m+n>2 ^ '' ' 

where 72m,n are appropriate derivatives of [3fh- We consider a truncated form of [3ujint{<t', ??), 
because otherwise analytical results for the fluctuation-contribution to the grand potential 
are not possible. Strictly speaking, the above expansion can be truncated for 0^0 and 
?7 — > 0. For given values of (f) and r/, in particular for the results of our calculations in the 
ordered phase, however, the truncated expansion may be oversimplified, especially for small 
values of pg and for large amplitudes of the fields. 

In the field theory the grand potential and the charge-density correlation function are 
given by 

fi = -A;TlogH, (14) 

and 

(0(x)0(xO) = S-i j D^j D0e"^^^""0(x)0(xO (15) 

respectively, where 

H= f D^ f Dct^e-^'^''^' . (16) 

n 

In the weighted-field approximation (WF) introduced in Ref . [13 , and described in more 

nnri 

detail in Refs. [2a,|2^l^, the field ^^(x) is approximated by its most probable form for each 
given field </'(x). Another words, for a given field 0(x), the field ^^(x) is determined by the 
minimum of /3Afi*^'^[0, 77] ((5/3Af2^''^^[0, r^J/tJr] = 0), and can be written in the form 

Wp(0(x)) = E^0(x)'", (17) 



3511. Insertion of 



where the coefficient a„ is given in terms of 72m j' such that m + j <n+l 
?7vyF(0(x)) into Eq.Q leads to simplified forms of Eqs. (fTT)|) and (|TH|l . 

D0e^^^^^'^[*l (18) 



and 

(</.(x)(/.(x')) = S-i / D0e"''^=^^['^V(x)0(x') (19) 

respectively, where 



oo 



/?^e//[0] = lff 0(x)C°,(x - x')0(x') + E ffe / <^'"(^)- (20) 

The coefficient A2m is given in terms of 72fc,n such that k + n <m J35|. For the fluctuation- 
contribution to the average density we obtain 

(r^(x)) = S-i /'D0r^H.F(0(x))e-^^^//[<^] = ^ ^(0(x)'"). (21) 

n 

Note that when the expansion in Eq.()20|) is truncated at the term oc 0^™, then in a consistent 
approximation the expansion in Eq. ()21|) should be truncated at the term oc {(f^^) with 
n < m — 1. Otherwise a„ would contain the coefficients 72^,^ that in (37ieff[4>] ^^^e not 
included. 

We should mention that on the sc lattice the WF approximation yields quite good results 
for the locations of the continuous order-disorder transition, and of the tcp 35[. In gen- 
eral, for the approximate WF theory we cannot expect exact dependence of the calculated 
quantities on Pq. 

In this work we shall limit ourselves to the 0^ theory, with Tie// approximated by TiwF 
of the form 



The explicit forms of A4, Aq in the WF theory are given in Appendix A for the CS reference 
system. The line of instability of 'Hwf[(P] (A-line) is given in Eq. (fTT|) . For stability reasons 
the expansion in Eq. ()20|) can be truncated at the term ~ 0^"", if the corresponding couplig 
constant is A2no > 0. In the case of the CS reference system we have ^4 > for pg > 
pI^ ^ 0.09795, and ^4 < for p^ < pl^. We find ^g > outside the density interval 
Ptc < Po < 0.1541. Negative coupling constants were also found for the RPM in Ref. J51|. 
We shall not calculate any quantities for pl^ < p^ < 0.1541, where the functional (j22p is 
unstable. Near the above range of densities our results are particularly strongly influenced 
by the lack of the terms 0(0^), and are less accurate than elsewhere. 



Note that C^Ak) given in Eq.® assumes a minimum for k = kb > 0, and can be written 
in the form 

C'^4k) = P*{To + AV{k)) (23) 

where 

f3'To = \ + /3*Vih) (24) 

and 

AV{k) = V{k) - V{kb) -k^k, V2{k - kbf + 0{{k - kb)'). (25) 

Near the hne of instabihty oiHwF, we have /?*ro -^ (see Eq. (|TT| ). Because C^^{k) assumes 
a minimum for k = kb, the fluctuations 0(k) with |k| ^ kb dominate, If the fluctuations with 
k significantly different from kb are irrelevant, i.e. for tq <^ V2kl [40], the term 0{{k — kb)') 
in Eq.(j2SI) can be neglected, and 

Cl^{k)^(3*To + (3*V2{k-kb)'. (26) 

I — I 

Eqs. fl26|) and (j^^ with Aq = are of a similar form as in the Brazovskii theory [40.]. In 

the next section we derive an approximate form for the grand potential in the 0^-theory 
(Eq.(j22I) and (PBj) ). by generalizing the Brazovskii approach. 

III. CONSTRUCTION OF THE BRAZOVSKII-TYPE APPROXIMATION FOR 
THE RPM 

In the charge-ordered phase characterized by a charge-density profile that is periodic in 
space, the fluctuating field can be written in the form 

0(x) = <|.(x) + ^(x), (27) 

where 

$(x) = (0(x)) (28) 

describes the ordered phase with a particular symmetry. In the ordered phase Eq. ()16|) can 
be rewritten in the form 

S = exp ( - PHeff m) [ Dtp exp ( - /37^/;„e[$, ^]) , (29) 



where 

HfiuA^, ^] = K//[$ + ^] - K//[$]. (30) 

For the (f)^ theory fEq. ()22|) ) we have 

^fU^, i^] = \l I ^(x)C,^;"^(x - x')^(x') + I Ci(x)^(x) + (31) 

*/ X •/ X •/ X 

i /c73(x)V^(x)3 + i /'c4(x)V^(x)^ + | /'c5(x)V^(x)^ + r^ /'^(x)^ 

t/X *t/X 't/X '</X 

where exphcit expressions for C^^"'^ and Cj are given in Appendix B. 

By inserting Eq. ()29|1 into Eq. ()14|1 we obtain a functional of the charge- density distribution 
$(x), 

(32) 



-/5fi[$(x)] = -/3K//[$(x)] + log [y Z^^exp (^ - /37^/,„e[$, M 

The above equation gives the grand potential in a system with the charge density constrained 
to have the form $(x). For the theory given by the coarse-grained Hamiltonian (PUJ) . or by 
its truncated version p2|l . this result is exact. For the charge distribution given by $(x), 
the first term in Eq. ()H2|1 can be directly calculated. In order to obtain an approximation for 
the fluctuation contribution, we rewrite l-ifiuc[^iip\ in the form 



Hfiuci^, ^] = 7^g[$, ^] + A7^[<l>, ij] 



(33) 



where Ticf"^,"^] is the Gaussian contribution, 

nci^M = iy"y"^^(x)C^^(x-x')^(x'). (34) 

(7^0 (x — x') is inverse to the exact charge-density correlation function, i.e. in Fourier repre- 
sentation C^^{k) = 1/G^^{k), where 



G^^{k) 



(35) 



Next we make an assumption that /3A7i[$, ip] can be treated as a small perturbation. When 
such an assumption is valid, we can write 



log 



/ 



Dtpe' 



-PCHg+^-H) 



log 



Dilje 



-mc 



i-l3/\n + o[{i3/\nf 



log / D^e-'^^^ + log 1 - {(3^n)G + 0{{l3l^n)l) 



(36) 
(37) 
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FIG. 2: Feynman diagrams contributing to C^^f, in the disordered phase, to two-loop order. Shaded 
circles and a bullet represent A^ and Aq respectively. Lines represent G^... In the self-consistent 
theory the lines represent G^i. 

where {■■■)g denotes averaging with the Bohzmann factor e"^^*^. Assuming again that 
{PAT-C)g is small, we obtain 



pn[^] = pHeffm - log f D^e-^'^^ + {pAn)G + OiipAH) 



(38) 



In the uniform phase jSQ is given by the same expression, but with $ = = /37ie//[0]. 
In practice the exact form of 0^,6 cannot be calculated analytically. In the perturbation 
, |53[ 6*00 is given by Feynman diagrams with the 2?T,-point vertices A2n- The 
vertices at x and x' are connected by lines representing G5,<a(x — x'), and all lines are 
paired. The corresponding expressions are integrated over all vertex points, or in Fourier 
representation over all G^Jk)-lme loops. In this work we shall follow the self consistent, 
one-loop Hartree approximation for C^^ |40|. The one-loop contribution to C*,^,^ (Fig.2a) is 
proportional to A4 J^ G^^{k). In the effectively one-loop cj)^ theory (j^ . another contribution 
to C^if,{k) is given by a diagram (Fig.2c) that is proportional to -4.6(/,^G°^(/c))^ 
The symmetry factors of the graphs are calculated according to standared rules 



52 



m 



35|. 



In 



the sel 



26 



35 



consistent, effectively one-loop approximation C^^{k) assumes the approximate form 

C^^{k) = r + f^*AV{k), 
where r = C^^{kb), and by using (jH^ . (j^ and (f7m) - (f7Sjl we obtain 



(39) 






V A\ 



V 



(40) 
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In the above V = J 1 is a volume of the system and 

G{r)^{i;{^f) = j^G%{k). (41) 

The remaining diagrams (including the one shown in Fig.2b) are negligible in the 0^ theory 
for Ai\/ f3*V2kb <^ r [40]. When the above condition is not satisified, the neglected diagrams, 
apart from a modification of the form of r, yield aditional, fc-dependent contribution to C¥j^ 
in Eq.(jnni)- Inclusion of such contributions goes beyond the scope of this work. 

In general, the integral in Eq. (j4H) cannot be calculated analytically. In fact the integral 
diverges because of the integrand behavior for k -^ oo. However, the contribution from 
A; — i> oo is unphysical (overlaping hard cores). When the fluctuations with k ^ kb dominate 
(r <^ l3*V2kl), then the main physical contribution to G{r) comes from k ^ kb- In this case 
the regularized integral is J40| 

/■ 1 r 1 _ 2aVT* 

^^^^ ~ k r + P*AV{k) "^''^^ L r + P*V2{k-kbr ~ ~7^' ^ ^ 

where 

a = kl/{47rV^,). (43) 

The Eqs.din} and (jl^ are to be solved selfconsistently for the ordered and the disordered 
phases. In the disordered phase, i.e. for $ = 0, r is denoted by rg. 

The second term in Eq. (jnH|) . with C^ff,{k) approximated by C^^{k) (see (|H^). is 

log / DV^e"^^^ ^ -2aVf*^V, (44) 

where the approximation ()25p and the same regularization as in the case of Eq. ()42|) were 
used. For the last term in Eq.ljHHD we find (see ((221), (ED), (123) and (gOD) 

^,^n)Mv . _£|)!^, _ £!!_)! ^. _ £M!^. I m. (45, 

The above results and Eq. ()42|l give the explicit form of the grand potential 



pn{pl r ; $; r]/V = PHwf{pI, T*; <^]/V + 2aV?F (46) 



3 \ r J A r Jy. V 

where r = r(pQ, T*; $] is a function of Po^ T* and a functional of $(x) that is to be determined 
from Eqs. (flUj) and ((I^. For given values of pg and T*, the value of the (dimensionless) 

12 



grand potential for a considered phase corresponds to the minimum of Pfi{pQ, T*;^;r] with 
respect to $(x), with pg and T* fixed. 

T* represents temperature, but pg differs from the average number- density when the 
fluctuations are included (see ()21|l ). The lowest-order fluctuation-induced density shift, given 
by the first term in Eq. (pT|) . yields the leading contribution to the average local density p*(x) 
of the form 

Pt(x) = pS-^(^(r) + $(x)2). (47) 

The above gives the density shift in the 0^ theory. Higher order terms in Eq.(j2^ can be 
included simultaneously with higher-order terms in Eq. ()22|l . In the 0^ theory the next-to- 
leading order term in Eq.(j2^ leads to the following approximation for the average density 



p;(x) = pt(x) + | 



<l>(x)^ + 6$(x)2e;(r) + 36;(r)2j, (48) 

where 02 is expressed in terms of 72m,n in Appendix A. The thermodynamic density is given 
by the space-averaged density profile according to 

where the integration Jy is over the unit cell of the ordered structure, and V^ is the unit-cell 
volume. Explicit expression for the average density in the liquid is given in Appendix A. 

In practice a determination of the equilibrium charge-density profile $(x) and the phase 
transition between the charge-ordered and charge-disordered phases from Eqs. (j46j) . (J4U|) and 
(j42p is difficult. The problem simplifies greatly when a form of $(x) is limited to a particular 
function that depends on several parameters. In this case the functionals are reduced to 
functions of several variables, and the problem of obtaining a minimum of f3Q for a given 
class of functions becomes tractable. 

For an ordered phase of a particular symmetry, $(x) can be written as a linear combina- 
tion of functions gi{x) forming a corresponding orthonormal basis J48|, 

$(x) = 5^<|.,(7,(x), (50) 

i 

where $4 are the corresponding amplitudes. In Fourier representation the basis functions 
can be written in the form 

^.(k) = 7= E (^(k - k^,) + 6{k + k^j) , (51) 
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where for the considered symmetry 2nj and k^^ are the number of vectors and the j-th vector 
in the i-th shell respectively. In order to specify the structure we need to know the vectors 
k]^ that determine the size of the unit cell of the structure, apart from the amplitudes $,. 
We assume that the vectors forming the first shell correspond to the wave-vectors of the 
most probable excitations in the uniform phase. In the theory outlined above such wave- 
vectors are determined by a minimum of C^^{k), since it yields a maximum of the probability 
oc exp(— /^TYg)- In the one-loop approximation C^^{k) assumes a minimum for |k| = kh, we 
thus assume m^^l = kf,. 

When the form of $(x) is restricted to Eq. ()5()|l and |k-J^| = ki,, then for given pg and T*, 
r = r(pQ, T*; $] and l3Q(pn, T*; $; r] become functions of the amplitudes $». Typically, only 
the first one- J40|, |47| or two shells J48| in Ea. (l5UD a re taken into account in studies of the 
fluctuation-induced first-order phase transitions 40|,|47|. Our explicit results are obtained 
in the one-shell approximation, 

$(x) = <l>i(7i(x). (52) 

For a few points we also considered the two-shell approximation, but the long formulas will 
not be given here. In the one-shell approximation we can write 

<I>(X)2 ^2 f H^y ^4 /•$(X)6 . 



where 

are the geometric factors associated with a particular symmetry of the ordered phase. In 
the one-shell approximation the above relations should be inserted into Eqs. ()40|) and 
and the extremum condition for f3Q can be written in the explicit form 



r+(^A, + aA^){^-l)^l + ^{..-5s.)<l>t = 0. 



(55) 



The resulting set of equations can be solved for each point (pg, T*) with respect to r and $i. 
When the results are inserted in Eq.(|46|). the minimum of the grand-potential is obtained 
for a given pair of So and k^, i.e. for a chosen structure of the ordered phase. The above 
method of obtaining the grand potential is equivalent to the method used in Ref . J40, |42| 
and outlined in Appendix D. We verified our calculations by comparing the results obtained 
by both methods. 
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IV. TRANSITION BETWEEN THE CHARGE-ORDERED AND CHARGE- 
DISORDERED PHASES 

Let us focus on effects of fluctuations on tlie A-line, which on the MF level is given in 
Eq. ljllj) . At the A-line the system becomes unstable with respect to the dominant charge- 
density fluctuations with the wave number kf, > 0. At the boundary of stability, the second 
functional-derivative of the grand-potential functional (|HSj) is C(i,^{kb) = 0. In the effectively 
one-loop self-consistent Hartree approximaiton we have C¥Jkh) = tq, where tq is a self- 
consitent solution of Eqs. ()4()|l and 1)411] with $ = 0. For tq ^ 1 the approximation (jl^ is 
valid, and we easily flnd that ro = is a solution of Eqs. ()4(jp and ()42j) with $ = only for 
T* = 0. Thus, in a presence of fluctuations the A-line disappears. 

For a given thermodynamic state there may exist one or several minima of f3Q (Ec\. ^^ ). 
associated with ordered phases with different symmetries. The lowest value of the grand 
potential for given pg, T* corresponds to the stable phase. Other phases are metastable (un- 
stable) if a minimum of the grand potential exists (does not exist). At the phase coexistence 
two minima of the thermodynamic potential, corresponding to the two coexisting phases, 
are equal; other minima, if present, are associated with larger values of the grand potential. 

When the expansion for $(x) in Eq. ()50|) is truncated at the flrst shell, analytical results 
for the phase coexistence are possible. However, only a very crude approximation for the 
ordered structure can be obtained. In this study we shall limit ourselves to analytical 
calculations in the one-shell approximation. More accurate results for the structure of the 
ordered phase can be obtained numerically in a future work. 

We are interested in a stability of the ionic crystal. In the case of the CsCl symmetry 
{PnSm), the flrst shell is formed by the three vectors k^/A^b = (1, 0, 0), ^l/kb = (0, 1, 0) and 
kfe/fcfe = (0, 0, 1), i.e. rii = 3, and in real space 

g[{x, y,z) = —^( cos{kbx) + cos{kby) + cos{hz)) . (56) 

The above flrst shell determines the so called P structure. Unfortunatelly, topological prop- 
erties of P differ from that of the CsCl crystal. Namely, the P structure is bicontinuous, 
i.e. the surface 5'i(x) = separates space into the positively and negatively charged regions, 
and both regions are continuous, as shown in Fig.3. In the ionic crystal, however, the posi- 
tively and negatively charged regions are topologically eqivalent to spheres separated by the 
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FIG. 3: The surface (7i (x) = in the unit cell of the periodic structure P. This surface separates 
the positively and negatively charged regions. The lattice constant is 211 /kh « 2.55 in ion-diameter 
units. 

uncharged solvent of nonvanishing volume. In addition, the nn distance in the P structure 
is v^7r/A;f, ^ 2.2. This distance is much larger than in the actual CsCl crystal. The nn 
distance in the ionic crystal is closer to the nn distance, vr/fcf, k, 1.27, in the case of the one 
dimensional ordering (lamellar phase), where ni = 1 and 

g[^'^{x) = V2cos{kbx). (57) 

Since the precise structure of the crystalline phase cannot be determined within the one- or 
two-shell approximation, we consider both phases to find and compare the transition lines 
between them and the disordered phase. In this way we gain some insight in the approximate 
location of the actual phase transition. 

A. MF approximation 

In the MF the fluctuation contribution to (J38|) is neglected. In the one-shell approximation 
the grand potential is a function of the amplitude $i, 

PUwFipo, T\ ^{)/V = ^/3Vo$? + ^so^t + ^^o^t, (58) 
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where the geometric factors So and Ko are found to be 



, 3/2 lamellar structure 
So={ ' (59) 

5/2 P structure , 



and 



, 5/2 lamellar structure 

no={ ' (60) 

155/18 P structure . 

For Pq > pI^ the order- disorder transition is continuous and coincides with the line of insta- 
bility PHI , whereas for p^ < p^^ the order-disorder transition is first order, and occurs when 
the condition 



dpnwF{po,T*,<!> 



- = = /3Hh^^(Po,T*,$i) (61) 



(9$! 

is satisfied. The expressions for the transition lines T^*^ and Tp, and the amplitude $i are 
given in Appendix C fEqs. ()75|) and ()76|) ). It turns out that the P phase is only metastable. 
However, the relative difference (T^*^ — 7p)/T;*^ is very small. 

The density in the ordered phase can be obtained from Eq. ()47|) by neglecting the fluctu- 
ation contribution. In the lowest nontrivial order we have 

P*{^)=P*o-^H^r. (62) 

270,2 

We verified that Eq. (jU^ yields p±(x) = (p*(x) ± $(x))/2 > for all space positions. The 
space-averaged density in the ordered phases is given in Eq. ^^ . The explicit expression for 
Ap* = p* — Po is given in Appendix C. The resulting density-temperature phase diagram is 
shown in Fig. 4. 

The diagram shown in Fig.4 is obtained with the functional j3AQ^'^^[(f), rj], Eq.©, approx- 
imated by the functional Hwf, Eq.(j221)- In order to verify the accuracy of the approximation 
(1^ . we calculated the functional (3An'^^^[^{x), p*{x) - p*] along the line T,;^(p*), for the 
fields $(a;) and p*{x) — pg that yield TiwF = 0. The result shown in Fig. 5 indicates that 
our approximate functional ()22j) yields the poorest accuracy in this part of the phase dia- 
gram where Aq is very small (see the discussion below Eqs. p3p and ((221))) and the term oc 0^ 
should be included. We also considered the two-shell approximation for $(x). We found very 
similar results, with somewhat lower transition temperatures, and with a smaller difference 
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FIG. 4: Density-temperature MF phase diagram obtained from the approximate functional (|58|) . 
Temparature T* and densityp* are in standard reduced units defined in sec. II. 

between them. We conclude that the transition-temperature obtained in the approximate 
theory (Eq.(j22I)) is overestimated. 

Our main concern in this work is to determine the fluctuation contribution to the grand 
potential in Eq.()38|). We shall not attempt to find better MF results by numerical mini- 
mization of the functional (|7j). 



B. Effects of fluctuations on phase transitions 



In this subsection we include the fluctuation contribution to f2 in Eq.( 
We consider two cases, the ^^ theory {Aio = in the above equations), as in the Brazovskii 
work |40|, and the 0^ theory. 



1. Results of the (j)'^ theory 

The (f)'^ theory is stable for p^ > p*^, and for such densities the term 0(0^) can be neglected. 
In this case analytical solutions for tq and r of Eqs. pnj) and (jKK|l can be obtained. Physical 
solution for r corresponds to the lowest value of $i. The transition lines between the uniform 
and the two ordered phases are shown in Fig. 6. In the one-shell approximation the lamellar 
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FIG. 5: The MF grand potential /5AQ^^[$(x), p*(x) - p^] (see Eq.©) along the approximate 
transition line (f7S)) . where TCwF = 0. $(x) and /0*(x) are given in Eq. (|5^ with ((THJ, and in Eq. 
(|H^ respectively. 

order turns out to be more stable than the P phase. We also considered the much more 
tedious two-shell approximation for a few points. We found $2 significantly smaller than 
$1, and the phase-transition lines shifted to somewhat lower temperatures compared to the 
one-shell approximation. The rest of our results is obtained in the much simpler one-shell 
approximation. 

Recall that our results rely on the approximate Eq. p^ . which is valid provided that the 
condition r -C f2fc^/^*(Po) is satisfied. We verified that along the coexistence lines T*(pq) 
(Fig. 6), the r and tq are one and two orders of magnitude smaller than V2kl/T*{p'^ k, 
6/T*(pq), respectively. Thus, close to the phase coexistence the approximation (jl^ used in 
our calculations is valid. However, the condition A4\/j3*V2kf, -C r, under which the disre- 
garded diagrams (including Fig. 2b) can be neglected H40| is satisfied only for high densities. 
Therefore the accuracy of our results increases with increasing density. 

In Fig. 7 the temperature at the transition to the lamellar phase is shown as a function 
of the most probable density (MF result) and as a function of the average density, given by 
the approximate expressions (jTfj) and (jEHI)- The average density is a nonmonotonic function 
of T* along the phase coexistence (Fig. 7) for T* < 0.15. As already discussed in sec. 3a, for 
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0.20 




FIG. 6: The fluctuation-induced first-order transition lines between the fused salt and the ordered 
phases in the (j)'^ theory in the (pg, T*) phase diagram. The transition to the P phase is metastable. 

the corresponding range of pg the approximate functional ()22j) is oversimplified. Moreover, 
the neglected diagrams (Fig. 2b) may yield a relevant contribution to the grand potential for 
low densities. 

Let us compare the temperatures at the continuous transition in MF and at the first-order 
transition in our theory (Figs. 4 and 7). In particular, for p* = 0.8 we find T* ^ 1.29 and 
T* ^ 0.43 in the first and in the second case respectively. For T* ^ 0.43, on the other hand, 
we find in MF the transition-density pg ^ 0.28, a much lower value than in our theory. As 
we see, in this approximation the fluctuation-induced shift of the liquid-phase boundary is 
substantial. However, for p* = 0.8 the temperature at the transition between liquid and the 

n 

CsCl crystal obtained in simulations |^ is T* ^ 0.1. Before we identify the charge-ordered 
phase, we need to find out how the transition temperature and density change when better 
approximations for TiwFi for the function $(x) and for the average density are made within 
our theory. In the next section we study the role of the (j)^ term in Eq. ( 
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FIG. 7: Temperature at the liquid- lamellar phase transition in the (/)^ theory as a function of the 
average density in the liquid phase at different levels of approximation in Eq. (|2H) . Thin solid line 
corresponds to the zeroth-order term in Ea. (|21() . i.e. p* is approximated by p^. Thick solid line is 
obtained by including the leading-order contribution to the fluctuation-induced density shift, i.e. 
p* is approximated by p\ fEa.H47|)). Explicit expression for the average density in the liquid is 
given Ea. (|M|) . 

2. Results of the (j)^ theory 

In this subsection we determine the effect of the 0^ term on the phase behavior. Analytical 
solutions for the phase transitions can be obtained by using the original Brazovskii method 



40L |47| , if in equations determining the phase transition the terms of the highest order in 



$1 are neglected. This is justified when $1^1. In Appendix D we explain the key steps of 
the calculations. The full equations in the 0^-theory can only be solved numerically. Results 
for the transition lines between the uniform and the two ordered phases are shown in Fig.8, 
where analytical results of the approximate theory and numerical results of the full 0^-theory 
are shown as lines and as symbols respectively. The liquid-phase boundary, T*(pq), is shown 
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FIG. 8: The fluctuation-induced first-order transition lines between the liquid- and the ordered 
phases in the ^^ theory in the (pg,T*) phase diagram. Lines are the analytical solutions of the 
approximate equations (Appendix D) and symbols are the numerical solutions of the full equations 
described in the text. The transition to the P phase is metastable. 



in Fig.9 as a function of the average density at different levels of approximation in Eq. (PT|) . 
The thick solid line is obtained from Eq.(j25 with the two leading-order terms included, 
i.e. in an approximation consistent with Eq.(j22I)- We verified that Eq . (j42j) used in our 
calculations is also valid in the 0^-theory. 

By comparing Figs. 6 and 8 we see that in the 0^-theory both transition lines are signif- 
icantly shifted to lower temperatures compared to the 0^-theory. Fig.9 shows that in the 
consistent approximation for the average density, higher density at the phase coexistence 
is obtained. The results of the 0^ theory are thus closer to the simulation results for the 
liquid-CsCl transition, but the transition temperatures are still too high. For example, for 
p* = 0.8 we have T* ^ 0.43 and T* ^ 0.28 in the 0^ and in the 0^ theory respectively, 
whereas for the liquid-CsCI transition T* ^ 0.1 for the same density y. Note that the MF 
results (see the discussion at the end of sec. 4a) suggest that the truncated functional 
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FIG. 9: Temperature at the first-order phase transition between liquid and the charge-ordered phase 
in the (j)^ theory, as a function of the average density ()21() at different levels of approximation. Thin 
solid line corresponds to p* approximated by the MF result, p^. Along the dashed line p* is given 
by the space- averaged leading-order fluctuation-contribution to the average density p\ ( Eg . (14711 ). 
Thick solid line corresponds to P2 ( Ea. H48() ). where the next-to- leading orderd term in (|21() is taken 
into account. 



leads to overestimated transition temperatures compared to the original functional ((Tj). We 
can expect that by addition of the term ~ cf)^ in Eq. ()22|) and the third-order term in Eq. (PT|l . 
lower temperatures and higher densities at the transition should be obtained. Better forms 
of the charge- density profile $(x) should lead to lower transition-temperatures as well, as 
indicated by the (partial) results that we obtained in the two-shell approximation. Thus, 
systematic improvement of the approximations that we make in explicit calculations, should 
lead to systematically decreasing transition temperature and increasing density, and our 
results should systematically approach the simulation results for the liquid - ionic-crystal 
transition J4|. 

For the two quite different charge-ordered phases the transitions to the disordered phase 
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FIG. 10: The density profiles of cations, p%{x), in the charge-ordered phase at the coexistence with 
the hquid phase for two different densities in the (j)^ theory. Dashed hne corresponds to pg = 0.5 
(p* ~ 0.6) and the solid line corresponds to pg = 0.7 (p* ~ 0.75). x is in cr-units and p^ is 
dimensionless. 

are not far from each other (Fig.8). It is thus plausible that the transitions to the other 
ordered structures, including the stable one, are located in the same part of the phase 
diagram. 

Let us focus on the density difference between the coexisitng phases. In the high- 
temperature part of the phase coexistence the density difference Ap* between the coexisting 
phases is rather small (see Ref.|j| and Fig.l). By using Eqs. ()49|) and (jlSj) we find a vanish- 
ingly small density difference between the coexisting liquid- and lamellar phases. This result 
probably follows from the poor, one-shell ansatz for the charge-density profile. The density 
profile of cations, p+(x) = (p2(^) + '^(3^))/2, where P2{x) is given in Eq. ipHj) . is shown in 
Fig. 10 for two densities at the coexistence with the liquid phase. Beyond the effectively 
one-loop approximation we expect some (probably weak) dependence of the unit cell of the 
crystal on density, but its determination requires further studies. The shape of the density 
profile in the charge-ordered phase resembles an average density profile in a crystal (but only 
in one direction). Of course in the one-shell approximation the crystalline structure cannot 
be reproduced accurately. Our results for the liquid phase are more accurate, because they 
do not depend on the form of $(x), which is the weakest point of our explicit results. 
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V. SUMMARY 

In this work we considered the fluctuation contribution to the grand-thermodynamic 
potential in Eq. (j32|) within the field-theoretic description of the RPM. Our main purpose 
was a determination of an order, location and nature of the transition between the charge- 
disordered and charge-ordered phases. We obtained an approximate expression fEgs- PUj) . 
PUJI and (|121)) for the grand-potential functional of the charge-density profile $(x) in the or- 
dered phase. Our approximation is based on the self-consistent, effectively one-loop Hartree 
approximation, applied to the 0^ theory that was derived for the RPM in Ref. 35|. 

We found that in the continuum-space RPM the A-line of continuous transitions disap- 
pears when the charge- density fluctuations are included. Instead, a first-order transition to 
a charge-ordered phase appears. The range of temperatures and densities at the first-order 
transition to the charge-ordered phase is similar to the range of temperatures and densities 
at the liquid - ionic-crystal transition found in simulations j^. In the charge-ordered phase 
the average charge- density exhibits oscillations with a period ~ 2.5o" (beyond the one-loop 
approximation the period may be slightly different), and with an amplitude ~ 0.5. Thus, 
our results strongly indicate that the fluctuation-induced first-order order- disorder transi- 
tion should be identified with ionic crystallization. We can conclude that the A-line found 
in different theories of the MF type is in fact a MF indication of ionic crystallization. 

The explicit results for the first-order phase transition were obtained analytically within 
the simplest, one-shell approximation for $(x) (see (jH^ ). for two ordered structures: one- 
dimensional, lamellar phase, and three-dimansional, P phase shown in Fig. 3. We do not 
expect that the structure of the charge-ordered phase can be correctly reproduced on this 
simple level of approximation. Our approximate result for the crystalline structure is only 
a very crude approximation. The transition lines for the two different structures, however, 
are located close to each other on the {p*,T*) phase diagram. It is thus very plausible 
that the transition line to the stable phase should be located near the two transition lines, 
and conclusions concerning the approximate location of the order- disorder transition are 
justified. 

In the effectively one-loop approximation we find the same lattice constant of the CsCl 
crystal 2TT/kfy for a range of densities. Beyond the one-loop approximation we expect a 
weak dependence of the lattice constant on density. This result and Fig. 10 suggest that the 
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number of defects (vacancies) in the crystal coexisitng with the hquid (fused salt) increases 
with decreasing T*. The crystal melts either when at relatively low T* many defects are 
present (low p*), or when there are only few defects, but T* is high. 

Our results show a reasonable agreement with simulations even for the very crude ap- 
proximation for $(x), and we verified that by addition of further terms to Eqs.(|22p. (|47p 
and (j^ a better agreement should be obtained. The accuracy of the results can be signif- 
icantly improved within the approach developed in this work by choosing a better ansatz 
for the form of $(x). Note that the fluctuation contribution to f3Q depends only on global 
characteristics of $(x), i.e. on integrals of "^(x)^", where n < 3 in the case of the 0^ 
theory. Only the MF contribution depends on a detailed shape of $(x) through the term 
J^J^, <l>(x)C^^(x - x')$(x') = /^<l(k)C'^^(k)l>(k'). Numerical determination of the struc- 
ture of the charge-ordered phase is thus possible within our apporach, but it goes beyond 
the scope of the present work. We conclude that the field-theoretic approach developed in 
this work is suitable for a description of ionic crystallization on a semiquantitative level. 
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VI. APPENDICES 

A. Coefficients A^^Aq and a2 in the WF approximation 

The coupling constants in the WF approximation are given in terms of the coefficients 

72m,n |35|, 

A = 74,0 -3^^^^, (63) 

70,2 

. , J-72,i)(-74,i) , J-72,i)^(-7o,3) , J-72,2)(-72,i)^ .„.x 

A = 76,0 - 15 15 3 45 2 ' (64) 

70,2 7o,2 7o,2 

and in the CS approximation they assume the explicit forms 



1 - 20s + lOs^ - 45^ 
pf{l + 4s + As^-As^ + s^ 



. i- — ^UCl -r XU,5 — '±■5 -r ,5 /„-\ 

•^4 = - ,.3,. , .. , ... ... , ..^ (65) 



and 
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where 



A, = '-^^ (66) 



W{s) = 3 - 84s + 360s^ + 26Us^ + 1701s^ - 8736s^ 

+ 11240s^ - 8304s^ + 3861s^ - 1164s^ + 240s^° - 36s" + 3s^^ 



The coefficient 02 in Eqs. (PT|) and (jIH)) is 

-nO 

02 J- 4,1 



2 4170,2' 

where 



(67) 



po 6-72,2-72,1 3-73,0-72,1^ . . 

r4,i = 74,1 -2 • 68 

7o,2 7o,2 

In the hquid phase the exphcit form of the average density (jTfj) is 

a(l _ s)4 /pT 

Pi=Po + ^*(i + 4, + 4^2 _ 4^3 + s4) V ^' ^^^^ 

where Eq. ()42|) . Eq.()40|) with $ = 0, and the CS reference system have been used. The 
exphcit form of pg can be obtained in the same way with the help of Eqs. ()67|) . (jUHjl . 



B. Explicit forms of the functions CJ>^'^ and Cj in the functional ()31j> 

Ci'ri^ - X ) = Cj,(x - x') + (^$^(x) + ^$^(x))5(x - x'), (70) 

Ci(x) = I $(y)C7;,(y - x) + :^$3(x) + :^$^(x), (71) 



and 



C3(x) = A$(x) + 4^$(x)^ (72) 





C4(x) = A + ^$(x)^ (73) 

CrM) = <5(x). (74) 
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C. Explicit expressions for the phase transitions in the MF approximation 

The Une of the first-order transition and the amplitude of the charge-density wave are 
given by 



and 



- 2 •^i^o 






respectively. The space-averaged density shift has the form 

27o,2A«!,' 



AP- = t:."'':''L . (77) 



D. Explicit expression for the grand-potential difference 

After a substitution of Eqs.(27) and (28) into Eq.(22), the Brazovskii's equation of state 

[22], 

sn 

h 



can be written as 

h= (Ok,) + M + :^)$(_k,) (78) 

^l"^ + U^) L L L ^^^' ^^^ ^' ^ k'")<^(k')*(k")<i>(k'") 
+4j- [ [ [ I I 5(k6 + k' + k" + k'" + k"" + k"")n'^(k*)' 

^- Jk' Jk" Jk'" Jk"" Jk'"" j 

where C°^(k) and Q are given in Eqs. (26) and (43) respectively. For C^^{k) (see Eq.(27)) 
we obtain Eq.dHHI) with ^ . 

For |$| ^ 1 we truncate the expansions in $ in Eqs. ()78j) and ()40|1 at the 0($'^) and 
0($^) terms respectively. As a result, the correlation function is given in Eq. ()4()|l with the 
term 0($^) neglected. The truncated Eq. ()78j) . and Eq. ()4()j) give for h the result 

h = r|.(-k,) + (a + ^) j I Hk^)^k,) (i|.(k, - ki - k2) - ^<l(k,)) . (79) 

In the one-shell approximation the explicit form of the equilibrium condition h = is 

r+(A + ^)&n$? = (80) 
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where a = a\fT , with a given in Eq. ()43|) . and 

fen = ^ I I ^l(ki)^l(k2)^i(k, - ki - k2) - \. (81) 

For both the lamelar and P structures hn = — l/(4ni) . 

The difference between the thermodynamic potential in the ordered and the disordered 



phases can be obtained in the same way as in Ref. 



m 



and we find 



An = 1^$4 + Qr, (82) 



where 






l\^^{r + ^^'), (83) 



2 2r3/2 - 2^/3*To - A^a a 



^. = V-- -- - ^, (84) 



and 

^dv={^ + ^ + ^0-)dr: V = A,V~r + A,a. (85) 

After inserting ()84|1 and ()85|1 into ()83|1 we obtain an integral which can be calculated 
analytically (e.c. with the help of Mathematica). In order to obtain r and $i at the 
coexistence, we solve Eq. ()4U|) supplemented with the condition ()80p. As a result, we arrive 
at the explicit expressions for /Wt and r, which were used in a determination of the phase 
diagram shown in Fig.9, but are too cumbersome to be presented here. 
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